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Abstract. We investigate the entanglement dynamics of bipartite quantum system 
between two qubits with the dissipativc environment. We begin with the standard 
Markovian master equation in the Lindblad form and the initial state which is 
prepared in the extended Werner-like state: p%g{Q)- Wc examine the conditions for 
entanglement sudden death (ESD) and calculate the corresponding ESD time by the 
Wootters' concurrence. We observe that ESD is determined by the parameters like 
the mean occupation number of the environment A^, amount of initial entanglement a, 
and the purity r. For A'^ = 0, we get the analytical expression of both ESD condition 
and ESD time. For A" > we give a theoretical analysis that ESD always occurs, and 
simulate the concurrence as a function of jot and one of the parameters A^, a, and r. 



PACS numbers: 03.65.Ud, 03.65.Yz, 03.67.Mn, 05.40.Ca 



The entanglement dynamics of bipartite quantum system: Towards entanglement sudden death2 



1. Introduction 

Entanglement is responsible for the most counterintuitive aspects of quantum mechanics 
[1], and motivated many philosophical discussions in the early days of quantum physics 
[2]. Recently it has been regarded as a resource for quantum information processing [3]. 
In fact, entanglement is one of the key ingredient for quantum teleportation [H [5], E], 
quantum cryptography [7] and is believed to be the origin of the power of quantum 
computers, etc. However, a quantum system used in quantum information processing 
inevitably interacts with the surrounding environments (or the thermal reservoirs), 
which takes the pure state of the quantum system into a mixed state [8]. Thus, it 
is an important subject analyzing the entanglement decay induced by the unavoidable 
interaction of the interested systems with the environment [9]. In one-party quantum 
system, this process is called decoherence, and various methods have been proposed 
to reduce this unexpected effect [10], [11], [121 [IS]. In multiparty systems with non-local 
quantum correlations much interest has been arisen in the dynamics of entanglement. 
For example, entanglement sudden death (ESD), which means that it disappears at 
finite time, was discovered by Yu and Eberly [M] [15] . It differs remarkably from the 
single qubit coherent evolution. The interesting phenomenon has been experimentally 
observed for entanglement photon pairs [16] and atomic ensembles |17j . 

ESD puts a limitation on the time when entanglement must be exploited. The 
evolution of the entanglement and ESD have been analyzed and various interesting 
results obtained [l8l[T9l[20l[2ll[2a[23l[2a[25l[26l[27l[2g Ref. [23] investigated the 
time evolution of entanglement of various entangled states of a two-qubit system. For 
four different initial states they analyzed the entanglement sudden death conditions and 
time. From the simulation they got the conclusion that the ESD always exists except 
for the vacuum reservoir. However, it is still a question when the ESD occurs and what 
ESD time is in the general cases and universal initial states. The aim of this paper is to 
discuss the above problems for the standard Markovian master equation in the Lindblad 
form and the initial states, the extended Werner-like states: Pab(O)- 

The paper is organized as follows. We first introduce the Wootters' concurrence 
[29] and the extended Werner-like states: PAsi^)- Section HI we give a standard 
Markovian master equation [301 131]. The master equation is equivalent to a first order 
coupled differential equations when the initial state is the extended Werner like state 
Pa_b(0)- sections IV, we analyze the ESD conditions and ESD time for the initial 
state P*b(0). Conclusions and prospective views are given in Section V. 

2. Concurrence and initial states 

A useful measure of entanglement is the Wootters' concurrence [29]. For a bipartite 
system described by the density matrix p, the concurrence C(p) is 




(1) 
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where Ai, A2, A3, and A4 are the eigenvalues (with Ai the largest one) of the "spin-flipped" 
density operator ^, and 

(2) 



c 



y "y /r \"y ^ ^y 

where p* denotes the complex conjugate of p and ay is the usual Pauli matrix. C ranges 
in magnitude from for a disentanglement state to 1 for a maximally entanglement 
state. 

The general structure of an "X" density matrix [TH [T5] is as follows 
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(3) 



with X, y, z, w real positive and u, v complex quantities. Such states are general enough 
to include states such as the Werner states, the Bell states, et ai. A remarkable aspect 
of the "X" states is that the initial "X" structure is maintained during the Lindblad 
master equation evolution. This particular form of the density matrix allows us to 
analytically express the concurrence at time t a.s W-. 



C^{t) = 2max{0, |m| 



'xw, \v\ 



/yz}. 



(4) 



In the present work, we will analyze in detail the two-qubit entanglement dynamics 
in Markovian environment starting from the initial "X" states defined in Eq. (3). We 
will examine the exactly ESD time and entanglement evolution for two kind of special 
states, the extended Werner-like states. 



I AS, 



(5) 



where r the purity of the initial states, Iab the 4x4 identity matrix and 

$AB = (cos(a)|10) + sin(a)|01))AB, (6) 

with a measuring the amount of initial entanglement. The state p*^(0) has the following 
form 
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(7) 



Obviously, the state in Eq. (5) is reduced to the standard Werner state when 
a = 7r/4, and the Werner-like state become totally mixed state for r = 0, while the 
well-known Bell state for r = 1. Note that, the extended Werner-like state Pab{^) 
contain both separate state and entangled state. According to Peres' criterion |32j . 
when , , „i „ < r < 1 the extended Werner-like state would be entangled, otherwise 

l+2|sin2Q!| — — 05 

it would be separated. Fig.l shows the entangled and separated areas. 
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3. The master equation 



The standard two-qubit Markovian master equation is the foUowing Lindblad form 

2 



dp 
'dt 
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HuJn 



{2a-pa+-a+a-p-pa+a-}+—- }_^{2a+pa--a-a+p-pa-a+}, (8) 



where = l/(e^'s^ — 1), the mean occupation number of the environment osciUators, 
7o is the spontaneous emission rate. As introduced before, we assume the initial state 
to be a "X" state. When substituting (3) into the master equation (8) we obtain the 
following first-order coupled differential equations, 
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,(9) 



(10) 



u(t) = -(l + 2iV)7ou(t), 
v{t) = -(l + 2A^)7ot;(t). 

The solution of the previous master equation can be found by solving the system of 
differential equations. The reduced density matrix elements x(t),y(t), z(t),w(t),u(t) 
and v(t) are given 
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(11) 



(i) iV > 

x{t) = ci^-c2r2(t)-c47oA^r(t), 

yit) = C,+C2T\t)+C3T{t), 

z{t) = c, + C2T^{t)-csT{t)-caoT{t), 

wit) = Ci^-C2r2(t) + C470(iV + l)r(t), 

uit) = c,T{t), 

v{t) = ceT{t) 

(ii) N = 

x{t) = d{T^{t), 

yit) = d2r{t) + d,T{t)-d,T\t), 

z{t) = -dsr{t)-d,T^it), 

w{t) = di-d2T{t) + d^T\t), ^ ' 

uit) = d{T{t), 

v{t) = deT{t) 

where T{t) = e~^^+^^^'^°*, T{t) = e'"'"^, and the coefficients ci, C2, C3, C4, C5,C6 in Eq. 
(11), and di,d2,d3,d4^,d5, d^ in Eq. (12) are determined by the corresponding initial 
conditions. 



4. Entanglement dynamics with the initial conditions 
4.1. N >0 

We now analyze the entanglement dynamics. Starting from the initial states p^si^) 
Eq. (7), the coefficients of Eq.(ll) are determined as 

_ N{N+l) 

'^1 ~ (2Ar+l)2 ' 

_ N{N+l) _ 1-r 

^2 (2Ar+l)2 4 ' 

^3 = 2(2jv+i)2 + i(cos^(«) - sin^(a)), ^^^^ 

■^4 = ~7o(2Ar+l)2' 

C5 = r sin(a) cos(q;), 
C6 = 

Then the concurrence of p*^(t) is 

C{p%it)) = 2max{0, \u{t)\ - ^/x{t)w{t)}. (14) 
So the ESD appears when, 

\u{t)\ - ^x{t)w{t) < ^ u^{t) - x{t)w{t) < 0. (15) 



Obviously 

t-^+00, u^(t) - x(t)w(t) = — < 0, 



and 



(1 - rY 

t = 0, u'(t) - x(t)w(t) = sin^(a) cos^(a) - ^ ^ ^ . (16) 

16 
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Figure 2. (Color onlinc)(r = 1, a = 7r/4)Plot of concurrence of C*(t) vs "N" and 7ot. 

Thus, ESD always occurs if the initial state is entangled, which means 14.2 1 sin 2a\ — — 
according to Peres' criterion. 

Note that Ref. [23] found that when the mean thermal photon number is not 
zero, in the thermal reservoir the entanglement sudden death always happen from the 
simulation. However here we get the same result from the theoretic analysis. In Fig. 
2, 3, and 4 we simulate the concurrence as a function of 7ot and one of the parameters 
A^, a, and r, respectively. Fig. 2 is the concurrence as function of 7ot and A^, fixing 
the purity r = 1 and initial degree of entanglement a = j, the Bell-like states. It shows 
that ESD time is affected by A^. The smaller A^, the longer its ESD time. In Fig. 3, we 
plot the concurrence as function of 7ot and initial entanglement a when r = 1. 

When a = I or a = ^ the initial states are reduced to Bell-like states. The ESD time 
is sensibly affected by a. Fig. 4 is the concurrence as function of 7ot and purity r. As 
exemplified before, when < r < ]^^2|sL(2a)| initial state is separated, so we choose 
r from 1/3 to 1. It shows that the larger the purity r the longer the ESD time. 

Here we study the entanglement dynamics of bipartite quantum system in the global 
environment effect. However, Ref. [26] had shown that under the classical niose effect 
entanglement may experience a sudden death process even if the local coherence of one 
participating particle is well preserved and the other one decays to zero asymptotically. 
How do the local thermal reservoirs influence the entanglement dynamics? Up to my 
knowledge, the master equation needs to be reconstructed and the local environment 
effect embodied by the spectral density of the thermal reservoir J{lj,T). The one- 
body decoherence dynamics was studied in [121 132] under the local environment effect. 
Whether does the multipartite entanglement dynamics under the local thermal reservoirs 
hold under the classical noise effect proved by Yu etc al [26]? We will study it in our 
further work. 
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Figure 4. (Color online)(Q! = 7r/4, iV = 0.25)Plot of concurrence C*(t) vs "r" and 79^. 

4.2. N = Q 

The coefficients in Eq. (12) have the form 

di = 1, 

d2 = 1, 

4 = -^-rsin2(a), 

J _ l-r \^') 

"4 — — , 

d^ = r sin(a) cos(a), 
4 = 0. 

Thus 

u\t)-xit)wit)<0 
^ T^{t) - j^T{t) + ^ - ^ sm\a) cos(a) > 0. 
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So the ESD occurs when, 

1 < r < -l + Vl + 4sin-(2c.) .,9. 

l + 2|sin(2a)| 2sm'^{2a) ' ^ ^ 

and the ESD time is 

7ot* = ln(l - r) - ln[2(l - r + r'^ sin2(2a))]. (20) 

When 7ot G [7ot*,+cx3), the concurrence C*^ = 0. In Fig. 5 we plot the ESD area. 
If the initial state is pure, i.e., r = 1, we find that x{t) = 0. Then the concurrence 
is C{p'^^{t)) = 2\u{t)\ = I sin(2a)|e~'''"*, which implies that the nonexistence of ESD. 
After we submitted our paper, a closely related work appeared in [28]. Ikram et al 
extended their former result [23] to more general state for the zero temperature limit 
and provided some discussions on the thermal case. Comparing with us, there are 
three main differences. Firstly, although both papers studied the same system, a two 
2-level atom system, different initial states had been chosen. These two initial states are 
independent, and the results are complementary. Secondly, for pure 2-qubit entangled 
states in a thermal environment, Ref. [28] says that "We see in plots of figure 3 that 
entanglement sudden death always happens for non-zero average photon number in the 
two cavities and entanglement sudden death time depends upon the initial preparation 
of the entangled states." Here, a theoretical analysis is provided by Eqs. (14, 15, 16). 
Finally, in the vacuum reservoir we give the sufficient and necessary condition Eq.(19) 
for the entanglement sudden death, and answer the question when the ESD occurs and 
what ESD time is. Furthermore, the ESD area relies on the initial purity r and initial 
entanglement a, see Fig. 5. Here we give a more detail examination and analysis. 

The ESD condition and corresponding time for initial state PabI^) at = are 
summarized in the following table. 

TABLE I. ESD condition and time for initial state p%{0) at A^ = 



a 


(0, 2n) 


Condition 


1 / ^ / -l+-yi+4sin2(2a) 
l+2|sin(2a)| ^ ^ 2sin2(2a) 


Time 


7or = ln(l - r) - ln[2(l - ^/r + r^ sm'^{2a))] 



5. Conclusions 

In summary, we have analyzed the interesting phenomenon of entanglement sudden 
death determined by the dimensionless parameters hwo/kBT^a, and r. Sufficient 
conditions for ESD have been given for initial state P*b(0). We examine the conditions 
for ESD and calculate the corresponding ESD time by the Wootters' concurrence. 
We observe that ESD is determined by the parameters like A^, amount of initial 
entanglement a, and the purity r. For A^ = 0, we get the analytical expressions of the 
ESD condition and ESD time. For A^ > we give a theoretical analysis that ESD always 
occurs, and simulates the concurrence as a function of 7ot and one of the parameters 
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Figure 5. (Color onlinc)Plot the ESD area of the extended Werner-hke states p%g{t) 
with = for a e [0, 2tt]. The red hne is r = ^^^^■'^i/l^ , and the blue line is 

1 

r 



1+2| sin(2Q)| • 



A^, a, and r. In fact, these above results can also be obtain for other "X" initial states, 
for example p*^(0) = r\'$)ABAB{'^\ + ^lAB, where "^ab = (cos(a)|00) +sin(a)|ll))AB. 

We analyze the dynamic behavior of entanglement in the open quantum system But 
our ultimate aim is to control the entanglement such that it be a resource for practical 
realization. Despite of the noticeable progresses of entanglement, many fundamental 
difficulties still remain. One of the problem is ESD due to the interactions between 
system and environment; the other is that as the iV-particle increased, the entanglement 
becomes arbitrarily small, and therefore useless as a resource [3l] . Here we indicate that 
preparing some initial states can help prolong the ESD time. However, we think that 
it is a long way how to design some effective control field to make the A^-particle large 
enough and protect the entanglement that make it practically useful. 
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